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SIGNIFICANCE  AND  EXPLANATION 


Data  smoothing  is  probably  the  most  difficult  practical  problem 
confronting  Approximation  Theory,  in  part  because  it  is  hard  to  define  a 
satisfactory  concept  of  "smoothness"  mathematically.  This  report  adds  yet 
another  technique  to  the  available  arsenal  of  methods.  The  technique  is 
closely  modelled  on  the  physical  processess  used  when  ship  lines  are  'faired' 
or  drafted  with  the  aid  of  wooden  splines.  In  particular,  it  uses  something 
akin  to  weight  removal  and  something  akin  to  local  flattening  in  order  to 
achieve  its  goal  of  a  curve  with  few  inflection  points  yet  close  to  the  given 
data. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MKC,  and  not  with  the  author  of  this  report. 


DATA  SMOOTHING 


Guang-Chang  Dong 

Data  smoothing  is  an  important  problem  in  practice,  such  as  ship  lofting,  car  shell 
lines  design,  air-plane  lines  model  fairing,  mechanical  contact  line  smoothing  (for 
example,  the  smoothing  for  cam  data  design),  etc.  See,  e.g.,  [1]  for  a  mathematical 
discussion  of  this  problem. 

In  this  paper,  a  clear  concept  and  some  effective  methods  for  data  smoothing  are 
given;  this  is  a  part  of  author's  technical  reports  on  mathematical  ship  lofting  (2),  [3]. 


1.  The  concept  of  data  smoothing. 

Through  a  set  of  plane  data  (x^,y^)  (1  <  i  <  n)  with  end  conditions  draw  the  usual 
cubic  spline 

{  x—  X  )  (  X  “  X  ) 

y  =  Vmi,i+i(*-xi)  — —  t°i(2xi+rx_xi)+ci+i(xi+i+x_2xi) 


i+1  i 


(1) 


(Xi  <  x  <  x.+1,  i-1 ,2, . . , ,n) 


where 


m 


i ,  i+1 


i 

i 


yA/6  . 


We  call  c1,...,cn  the  coefficients  of  the  spline. 

We  think  of  the  spline  as  "smooth"  to  the  extent  that  the  curvature  of  the  spline  is 
varying  slowly  or  the  c^  ( i«1 , 2, . . . ,n)  are  varying  slowly.  We  make  this  more  precise  in 
both  a  qualitative  and  a  quantitative  way. 

Qu..jLitative  part.  One  requirement  for  a  "smooth"  spline  is  that  the  signature  of 
c^  should  change  slowly  or  the  number  of  inflection  points  R  of  the  spline  curve  should 
not  be  large.  If  the  original  R  is  large,  we  must  change  the  data  by  a  little  amount  so 
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as  to  reduce  R(  and  this  is  a  part  of  the  smoothing  procedure.  But  what  is  the  suitable 
final  value  for  R?  There  is  no  absolute  standard.  The  suitable  final  value  for  R  is 
dependent  on  our  requirement,  i.e.,  we  should  give  the  required  value  for  R  before 
smoothing.  Auxiliary  before  smoothing  we  also  give  the  qualitative  end  requirements:  He 
require  that  the  spline  on  left  end  is  convex  (given  D^l )  or  concave  (given 
The  same  for  right  end  (given  D2  =  1  or  D2«-1). 

Summing  up,  the  three  numbers  Dj,  D2,  R  are  given  at  the  outset. 

Quantitative  part.  Another  requirement  for  a  "smooth"  spline  derives  from  an  analysis 
of  the  force  from  a  duck  acting  on  a  wooden  spline.  The  exterior  force  on  node  (x^.y^) 
in  the  y  direction  equivalent  to  the  duck  holding  the  wooden  spline  fixed  is 
approximately  equal  to 


P 


i 


6EIei  := 


6EI  (' 


i 


i 


1-1 

i-1 


) 


(2) 


where  El  is  the  bending  rigidity  of  the  wooden  spline.  The  ideal  smooth  case  for  a  spline 
would  be  for  the  exterior  force  acting  on  consecutive  nodes  to  have  the  same  signature,  or 


Vl+1  >  0  (i-2,3 . n-2)  .  (3) 

Express  (3)  geometrically.  The  c(x)  :*  y"(x)/6  curve  consists  of  some  convex  and 
concave  polygons.  Then  (3)  demands  that  the  connecting  point  (x^,c^)  of  neighboring 
convex  and  concave  polygons  have  the  property  that  (xi-i<ci-1*r  *xi,ci>»  *xi+1'ci+l)  lie 
on  a  straight  line. 

We  do  not  attempt  to  take  the  ideal  case  as  the  goal  for  smoothing,  otherwise  one 
would  need  to  change  the  data  too  much.  Obviously,  keeping  change  of  data  to  a  minimum  is 
also  a  requirement  for  data  smoothing.  The  non  smooth  case  is 

ejej+1  <  0  <4> 


-2- 


with  both  e^  and  e^+1  large  in  absolute  value.  More  accurately,  the  amount  of  spring- 
back  when  we  release  a  single  duck  on  or  (xj+i<Yj+i>  are  large.  Hence  a 

part  of  our  smoothing  procedure  is  to  reduce  the  quantity  of  spring-back  to  a  certain 
amount  when  (4)  is  true. 

Before  turning  to  describe  the  methods  of  smoothing,  some  general  remarks  should  be 
given. 

1.  In  our  methods  of  smoothing  the  two  end  data  are  always  kept  fixed. 

2.  We  use  just  two  types  of  end  conditions.  For  the  left  end,  we  give  y^’  (fixed 
end)  or  else  give  no  slope  condition  (free  end).  I"  the  free  end  case,  we  always  add  the 
condition  c1  =  c2,  i.e.,  we  require  that  the  first  section  of  the  spline  is  a  parabolic 
segment.  The  right  end  is  treated  analogously. 

3.  In  addition  to  keeping  the  end  data  fixed,  we  can  fix  a  few  other  data  also. 

4.  If  the  spline  drawn  through  the  given  data  set  becomes  |y’(x)|  >  1  on  some 
segment,  then  the  principal  part  of  the  exterior  force  is  not  in  the  y  direction 
expressed  by  (2),  and  the  effect  of  our  method  comes  into  doubt.  For  this  reason,  we 
restrict  our  method  to  data  which  satisfy  the  condition  ly’(x)|  4  1  (x  4x4  x^). 

5.  We  can  treat  nearby  pairs  of  points;  i.e.,  pairs  of  points  (Xj,yj),  (x^+,,yj+1) 
with  xj+i-xj  less  than  a  quantity  given  at  our  disposal. 

2. _ The_Resilience_Method± 

Given  a  set  of  data  (x^.y^)  (1  4  i  4  n)  with  end  conditions,  draw  a  spline  with 

coefficients  c1t...,cn.  Change  y^  to  y^+Sy^  (2  4  j  4  n-1),  draw  another  spline  with 

coefficients  c. ,...,c  .  We  have 
1  n 


f 

I 


c, 

i 


rc^+c 


(j) 

i 


5yj 


(1  4  i  4  n) 


(5) 


where  c^^  <1  4  i  4  n) 
for  yj=1<  yj*0  (i+j) 


are  the  coefficients  of  the  j-th  unit  spline,  i.e., 
and  homogeneous  end  conditions.  We  can  prove  that 
(-1  ) ^  <0  (a  4  i  4  g) 


the  spline 

(6) 


I 

'  j 


'  I 
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where  a-1,  0»n  for  fixed  end  and  a-2,  S»n-1  for  free  end.  Moreover 


|c[j)|  diminishes  monotonely  and  quickly  (at  least  as  fast  as  2  ^ 

for  |i-j!  >2  . 

When  we  release  a  duck  acting  at  (Xj.yj)  on  a  wooden  spline,  i.e.,  we  set 


(7) 


=  66Ie ^  -  0 


then,  by  (5),  we  have 


5y  ,  -e  /e*^ 

Yj  y j 


(8) 


for  the  resulting  change  6y^ .  From  (6)  we  have  ej^  >  0,  so  that  by  (8)  we  have 
6y^e^  <  0.  This  is  the  source  of  the  name  "Resilience  Method”. 

Smoothing  ^rom  ejej+i  <  0  to  ejej+^  “  ®  by  only  releasing  one  duck  may  make  the 
datum  change  by  a  large  amount,  so  vie  consider  the  change  of  two  data  6y_,  and  6y_.  +  1>  We 
have 


c±  -  ci  +  cjj^ 5y j  +  cij+1  >5yj+1  f1  <  i  <  n) 


Let  Sy.  “  -6y.  .  and  take  the  minimum  value  of  |Sy.|  such  that  e.e.  ,  =  0.  It  is 
j  j+1  3  3  3+1 

easy  to  obtain  the  formula 


|e  , 


Ini! 


Syj  m  -6yj+i  ■  0 j , j+i  !  '  (’8ign  V  min{jj>_j3+v>'  ;(3+n_jj> 


}  . 


(9) 


ej  _ej 


ej+1  'ej+1 


In  the  actual  smoothing  process,  we  do  not  need  to  make  the  resilient  quantity 
disappear,  but  only  reduce  it.  We  distinguish  three  cases  as  follows: 
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1  •  when 


Sj-1Cj  >  °'  ej6j+1 


<  0. 


e3+16j+2  >  ° 


we  take 


-iy  . 

j+1 


0.8o  . 

3/3+1 


(10) 


2.  When  k  >  j+1  and 


6j-1ej  ” 


0, 


Vj+1  <  °'~ 


'Vl*k  <  °'  Ck' 


k+1 


>  0 


we  take 


6 y  *  0.4a  .  6y  *  0*4(-c  +c  )>•••* 

yj  3.  J+1'  ”  j+1  3,  j+1  j+1.  j+2  '  ' 

• 

6y,  .  *>  0.4(-o.  .  .  ,+o,  .  .) ,  4y,  =  -0.4o.  .  , 

k-1  k-2,k-1  k-1,k  Jk  k-1  ,k 


(11  ) 


3.  when  there  is  an  inflection  point  in  *xj'xj+l*'  and  in  *xj-i»xj+2*  there  is 
only  this  one  inflection  point,  i.e.,  when 


CjCj+1  <  °'  C3-1Cj  "  °'  Cj+1 C j+2  >  ° 


in  this  case  we  replace  a  by  0  in  (10)  or  (11). 

j ,  3+1 

II  we  meet  a  fixed  point,  then  at  this  point  we  simply  do  not  make  any  adjustment. 

This  process  is  called  one-time  relaxation.  How  many  times  shall  we  carry  out  our 
process?  When  we  find  that  m^xlo..  is  less  than  am  assigned  quantity  e  of  our 

choice,  we  carry  out  the  process  one  more  time  and  then  stop.  For  example,  for  ship 
lofting  we  always  take  e=3  mm. 


-5- 


Before  we  apply  the  Resilience  Method,  we  first  treat  all  'nearby'  points  as 


=  cj+1  =  S  mintlc  |,  |c^+2l)/5  (16) 

and  solve  for  Sy^,  ^yj+i  ■'■n  and  adjust.  We  get  the  final  result  (13). 

In  the  first  step  c_,  =  c^+1  =0,  i.e.,  y(x)  has  straight  line  segment  in 

[Xj,Xj+i).  This  corresponds  to  flattening  a  wooden  spline  by  pushing  a  ruler  against  it  so 
as  to  remove  the  Sc(x)  <  0  part  on  the  Bpline  curve.  This  is  the  source  of  the  name 
"Rule  Modelling  Method".  Our  final  result  is  c^  =  c_.  +  ^  expressed  by  (16),  i.e.,  y(x) 
has  parabolic  segment  in  *xj,xj+1^* 

We  now  give  a  detailed  description  of  the  Rule  Modelling  Method. 

Given  the  three  numbers  ,  D 2,  R.  Restrict  D ^  to  only  take  the  values  1,  -1, 
or  0.  For  the  left  end  we  require  that 

D  c  >  0  ( when  D  4  0 ) 

'  ( 17 ) 

C1°2  >  ®  (when  D1  =  0) 

similarly  for  the  right  end.  And  we  require  that  the  total  number  of  points  of  inflection 
of  the  spline  does  not  exceed  R. 

Suppose  the  above  three  requirements  are  not  all  satisfied.  We  start  by  finding  the 
interval  [Xj,xk)  for  adjustment  and  a  number  S  (S  =  1  or  -1)  satisfying 

Sc.  <  0,  Sc  <  0 ,  •  • . ,  Sc  <■  0,  Sc  >  0  (when  j  4  1),  Sc  >  0  (18) 

J  J'l  X  Kr 1 

(when  k  4  n)  . 

If  the  left  end  requirement  (17)  is  not  satisfied,  let  Xj  =  x^.  Let  S  *  D1  or 
S  =  sign  c2  according  to  4  0  or  else.  Determine  xk  by  the  conditions 
Sch  <  0  (1  <  h  <  k),  Sck+,  >  0. 

If  the  right  end  condition  is  not  satisfied,  find  S  and  [Xj,xkJ  similarly. 

If  end  conditions  are  satisfied,  and  condition  on  inflection  points  is  not,  then 
firstly  we  try  to  find  j  in  (2,  n  -  1]  such  that  c^_^c^  <  0,  and  Cj_.,Cj+1  >  0  and, 
if  such  j  exists,  let  S  =  sign  cj_^  and  xk  =  x^.  Otherwise  try  to  find  j  in 


-7- 


[2,  n  -  2]  such  that  cj-icj  *  ®»  *  0,  cj-icj+2  >  0*  and,  if  such  j  exists, 

let  S  -  sign  c^..,  and  xk  -  Xj+1,  etc. 

Thus  we  find  (Xj,xk)  and  S  satisfying  (18).  The  cases  k-j,  k-j+1,  k  >  j+1  are 
called  incorrect  small,  large,  and  extremely  large  wave. 

We  are  restricted  to  at  most  one  fixed  point  of  the  spline  in  the  interval  of 
adjustment  [x. ,xk) .  But  there  is  an  exceptional  case,  viz,  j=1,  the  signature  of  c  is 
not  correct,  we  need  to  adjust  and  y1 '  is  given.  In  this  case,  we  cannot  admit  any  fixed 
point  in  [Xj,xk]j  similar  for  the  right  end  case. 

Removing  an  incorrect  extremely  large  wave  can  be  reduced  to  the  case  of  removing  an 
incorrect  large  wave  by  deleting  the  nodes  between  x^,  xk  and  forming  a  new  spline  with 
n-k+j+1  nodes,  adjust  and  then  recover  the  spline  with  n  nodes  by  spline  interpolation. 
Thus  we  only  give  the  formulas  for  removing  incorrect  small  and  large  waves  in  the  middle 
or  near  the  left  end,  because  near  the  right  end  situation  is  similar. 

Formula  for  removing  incorrect  small  wave: 

when  1  <  j»k  <  n,  we  use  the  formula 


c.  -  v6yj-icji"1>+6yjCj3,+5yj+1c5j+1> 


Cj+o(ac‘  j'1)+Bc‘  j)nc‘: j+1)) 


(19) 


where  in  general  we  have 


1,  B  -  -1.5,  y  =  1 


(20) 


i.e.,  when  (Xj_1,y^1),  (^.y^),  (Xj+1  ,yj+1 )  are  not  fixed  points,  then  (20)  is  true. 
If  is  a  fixed  point,  let  a  «  0,  etc. 

First  step.  Let  c^  *  0  in  (19)  to  find  a,  then  adjust.  If  ScL  1  >  0, 

Sc^+1  >  0  are  true,  go  to  second  step,  otherwise  we  have  met  an  unusual  situation  and 
discuss  it  later. 
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From  Sc^_1  >  0,  Cj 


>  0  we  take 


Second  step. 


=  0, 


Scj+1 


c_.  =  Smin(|c  |,  ) c_.+1  )  J/5 

and  determine  a  by  (19),  then  adjust.  We  get  the  final  result 


(21) 


Sc(x)  >  0  (x.  _<x<x.  .)  . 

1-1  3+1 

When  j  =  k  =  1,  take  a=8=0,Y=1  in  (19).  Instead  of  (21),  we  take 


V8 


(22) 


Formulas  for  removing  incorrect  large  wave; 

This  is  the  case  (12).  If  there  is  no  fixed  point  in  ^*j'xj+l^'  adjust  as  in  (14), 
(15),  (16). 

If  one  of  (Xj,yj),  (x^+1,y^+1)  is  a  fixed  point,  consider  the  formula 


c . 

l 


ci+ci  6yj-1 


X  (j)*  ^  <j+1), 

+CJ  Oy.+C^  6y 


j+1^[j+2)6yj+2 


(1  <  i  4  n) 


(23) 


taking 


«yj_1  =  oa,  «y..+2  =  °8'  «yj+1  =  t5  •  (24) 

Take  y  =  1,  6  =  0  when  (Xj+,  ,y_j+., )  is  fixed.  Take  y  =  0,  6  =  1  when  (x^,yj)  is 
fixed.  Take  a  *  0  if  (x^_?,y^_j)  is  a  fixed  point  or  j  *  2,  or  Xj+1  is  an  extended 
end  point  (see  below)  which  we  marked  before,  otherwise  take  o  «  1.  Take  8=0  if 


r 


T 


<xj+2'yj+2*  is  a  fixed  point  or  j  -  n-2  or  Xj  ia  an  extended  end  point,  otherwise  take 
8-1. 

Solve  for  o,x  by  substituting  (23)  -  (24)  into  c  -  cj+1  "  i-hen  adjust.  This 
is  the  first  step.  Use  (16)  for  the  second  step  adjustment. 

If  there  is  a  fixed  point  (x*,y»)  satisfying  <  x*  <  x^+1,  and  k  >  j+1  and  we 

delete  the  data  points  between  (x^x^)  on  the  spline,  then  the  fixed  point  would  be  lost, 
too.  In  this  case  the  modeling  rule  cannot  move  freely  and  can  only  change  the  slope 
through  the  point  (x*,y*).  Following  this  idea,  we  give  the  adjusting  steps  as  follows: 

First  step.  Take  y*+y* (x*) (x^-x*) ,  y*+y * (x* ) (x^+^ )  instead  of  y^,  yj+1,  draw 
spline.  In  general,  this  spline  does  not  pass  through  the  point  (x*,y*).  Take 


y 


x^-x* 

xi+rxj 


, « 


Vi 

Vrxj 


in  (23),  take  a, (5  as  before.  Solve  for  o,T  by  substituting  (23)  -  (24)  in 
c^  «  c_j+1  *  0,  adjust,  draw  the  spline.  This  spline  contains  a  line  segment  in 
Xj  <  x  <  x  and  (x*,y*)  lies  on  this  line  segment. 

Second  step.  Denote  the  number  in  formula  (16)  by  X,  add  3X(x*-x^.)  (x^+1~x*)  to 
yj»  yj+1,  draw  a  spline.  The  final  spline  passes  through  the  point  (x*,y*)  and  contains 
a  parabolic  segment  with  y"  -  6X  in  x^  <  x  <  xj+i  and  satisfying  (13). 

If  j  -  1 ,  there  is  an  incorrect  large  wave  at  the  left  end.  If  there  is  no  fixed 
point  in  (x.j,x2)  and  y^'  is  not  fixed,  replace  (14)  by 


c^  -  Cj+c^1  *  6y^  *+cj2 'iy2  (1  <  i  <  n)  (25) 

where  c*1'  are  the  coefficients  of  the  unit  spline  given  by  y^  =  0  (1  <  i  <  n), 

V  - 


1,  and  homogeneous  right  end  condition.  Replace  (16)  by 


(26) 


c3/8 


Here  we  already  assume  that  SOj  >  0,  otherwise  we  shou) -  delete  nodes  to  treat  the 
incorrect  extremely  large  end  wave. 

Suppose  there  is  a  fixed  point  (x*,y*)  in  x1  <  x*  <  x^  or  y1 '  is  fixed.  In 
these  cases  the  modeling  rule  is  fixed,  it  is  the  straight  line  through  (x^y^),  (x*,y*) 
in  the  former  case  and  through  ( x1 ,y1 )  with  slope  y 1 1  in  the  latter  case.  Following 
this  idea  we  give  the  steps  for  adjustment. 

y*-y,  y'-y, 

First  step.  Replace  y^  and  y2  by  Xt^.x'  and  +  x«'_x  ^X2~X1  *  in  the  tomer 

and  replace  y2  by  y.j+y^ 1  ( x2~x1 )  in  the  latter  case.  Draw  the  spline,  in  general  this 

spline  does  not  pass  through  the  point  (x*»y*).  Replace  (14)  by 

ci  -  ci+ci3^y3  (1  <  i  <  n)  .  (27) 

Let  c^*0  in  (27),  adjust.  It  is  easy  to  prove  that  «  0.  Hence  the  spline 

contains  a  straight  line  segment  in  [x3,x2l  passing  through  the  point  (x*,y*). 

Second  step.  Add  31(x(-*‘),  3X(x2-x,|){x2-x*)  (where  \  -  c^/B)  to  y^*  and  y2, 
draw  spline,  Let  c 2  «  c3/8  in  (26),  adjust.  Then  we  can  prove  that  c^  -  c.j/8  and 
satisfying  Sc(x)  >0  in  x,)  <  x  S  Xj. 

Thus,  all  the  cases  for  adjustment  are  stated. 

The  rule  for  treating  an  unusual  situation  is  as  follows: 

When  we  do  the  first  step  for  adjustment,  if  at  least  one  of  Scj_i  4  Scj+h  *  ° 

(h  »  1  or  2  according  to  k  =  j  or  k  =  j+1 )  happens,  we  distinguish  four  cases. 

First  case.  Both  SCj  ,  <  0  and  Sc...  <0  are  true. 

-  j-1  1+h 

In  this  case  simply  turn  to  second  step,  only  need  a  changed  signature  in  the  right 
hand  side  of  (21)  or  (16). 
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Second  case.  One  of  Sc.  .  <  0,  Sc.,.  <  0  is  true,  with  c.  ,c.  _  >  0, 

-  3-1  3+h  ]-2 

°j+hCj+h+1  >  °* 

In  this  case  simply  turn  to  second  step  with  the  following  change.  In  (16),  replace 
S  by  sign  cq  when 


min{  lc  j_.j  I ,  lcj+1l)  “  lcal  <‘ 


similar  for  (21). 

Third  case.  One  of  Sc .  .  <  0 ,  Sc . , .  <  0  is  true,  with  Sc.  „  >  0,  Sc.  ,  ,>0. 
-  3-1  j+h  3-2  3+h+1 

In  this  case  we  need  to  extend  the  adjusting  interval.  After  a  new  first  step  no 

matter  what  the  signature  Sc.  ,,  Sc.  .  are,  we  never  extend  adjusting  interval  again 

3-1  3+h 

except  when  j  =  1  and  S  =  Dj.  Turn  to  second  step,  with  S  replaced  by  sign  c^, 

where  c  is  defined  by  (28). 
a 

Fourth  case.  When  Scj_1  *  0  and  **j'yj^'  ^xj-1'yj-1*  are  nearby  points  or 
Sc^+h  <  0  and  ( Xj+h_.j ,yj+h_1 ) ,  <xj+h,yj+h^  are  nearby  points,  in  any  of  these  cases  we 
always  extend  the  adjusting  interval. 

Remark.  In  the  above  cases,  when  j  =  1,  and  so  cj_1  does  not  exist,  when  j  =  2, 
and  so  cj_2  does  not  exist,  etc,  we  use  the  convention  that  then  the  corresponding 
requirements  Scj-i  <  0 ,  c^^c  2  <0,  ...  are  not  considered. 

It  only  happens  when  we  wish  to  remove  an  incorrect  end  wave  that  we  need  to  extend 
the  adjusting  interval  again  and  again  until  the  incorrect  wave  has  been  removed,  or  else 
we  arrive  at  the  extreme  case  [x^x^]  =  (x^,xn).  The  extreme  case  happens  when  one  of  the 
given  numbers  D^,  D^,  R  is  not  reasonable. 

In  principle  we  can  admit  in  (Xj.x^l  more  than  one  fixed  point,  in  this  case  the 
modeling  rule  is  fixed,  but  in  our  practice  we  have  not  met  it  and  we  do  not  want  to  make 


our  program  more  complex,  so  that  this  case  has  not  been  considered.  Only  a  similar  case 
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happened  in  practice  and  is  considered  by  us:  There  exists  an  incorrect  left  end  small 
wave  and  both  y, '  and  (x2#y2)  are  fixed,  and  in  this  case  we  adjust  by  simply  adding 
-l.lc^/c*^*  to  yj,  similar  for  the  right  end. 

We  have  finished  our  statement  of  removing  incorrect  wave.  If  the  three  requirements 
Dy,  Dj ,  R  are  still  not  satisfied,  then  we  adjust  cyclically. 

With  two  or  more  large  (and  extremely  large)  waves  existing  and  R  <  2,  we  cannot 
determine  from  D^,  D^,  R  which  one  is  the  incorrect  large  (or  extreme  large)  wave.  When 
this  situation  happens  (in  fact  it  happens  rarely),  we  print  c^( i-1 , . . . ,n)  and  let  a  man 
judge,  he  would  specify  the  interval  [ , x^l  of  incorrect  wave,  and  then  do  the  business 
of  removing  incorrect  wave  by  computer. 

Removing  incorrect  large  or  extreme  large  wave  may  produce  small  incorrect  wave. 
Removing  incorrect  internal  wave  may  produce  incorrect  end  wave.  But  the  probability  of 
this  happening  is  small  by  (7).  So  that  when  we  finish  the  process  for  removing  incorrect 
waves  once,  and  the  three  requirements  for  ,  D2,  R  are  still  not  satisfied,  we  should 
repeat  the  entire  process  again  and  again,  until  the  three  requirements  are  satisfied.  If 
after  repeating  the  whole  process  n  times  the  three  requirements  still  are  not  satisfied 
then  we  stop.  In  effect,  we  then  judge  the  three  given  numbers  d,,  D2»  R  to  be 
inconsistent  with  the  given  data. 

After  the  qualitative  requirement  on  the  spline  is  already  satisfied  and  should  we 
then  be  forbidden  to  use  the  whole  line  Resilient  Method,  in  this  situation  we  can  use  some 
supplementary  methods  of  local  adjustment  for  smoothing.  They  are  as  follows: 

A.  Adjustment  of  local  part  of  spline  which  is  too  flat. 

Assign  a  constant  K  greater  than  1,  usually  we  take 

K  £  {5,  3,  2.5,  2}  .  (29) 
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have  the  same  signature  and 


Case  1 . 


Cj'  °3+1 


K|c^|  <  min{ lc I ,  lcj+1li  •  (30) 

Case  2.  cj_i»  cj»  cj+i«  cj+2  heve  the  same  signature  and 

KlCjl  <  |c^_, i ,  K|cj+1|  <  |  cj+2|  .  (31) 

These  two  cases  are  called  locally  too  flat. 

For  case  1,  take  Sy^,  *yj+i  *nd  °»  S<  1  as  in  (19),  (20)  and  require  that 

K| 0^  |  «  mint  lc  I  .  lcj+1U  •  (32) 

We  can  find  a  and  adjust. 

For  case  2,  assume  that  there  can  be  at  most  one  fixed  point  between  ( xj-i '^j-l ) ' 
<Xj#y.j)  end  similarly  for  (*^+1  ,y.j+1  ) ,  <xj+2'yj+2)  *  Talte 

oy_)_l  -  oo,  Sy^  =  So,  6y..+  1  ■  yx,  5yj+2  “  (33) 

in  (19),  and  in  general  let 

a  =  6  =  2 ,  B  -  y  *  -1  .  (34) 


But  when  (Xj_1fyj_^)  is  a  fixed  point  let  a  »  0,  similar  for  B,  y,  5.  We  require  that 


Kc, 


'3-3' 


KC 


3+1 


'3+2 


(35) 


We  can  find  a,  x 


and  adjust. 


4 
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B.  Adjusting  of  local  part  of  spline  which  is  too  bent. 

Case  1.  (j-2<h<j+2)  have  the  same  signature  and 

I c j I  >  K  aax{|c  lcj+1l>  ,  (36) 


K  as  in  ( 29 )  . 

Case  2.  (j-2<h<j+3) 


have  the  same  signature  and 


|Cj |  >  t  |Cj+^|  >  K| Cj+2 I  • 


(37) 


These  two  cases  are  called  locally  too  bent.  Notice  that  in  some  cases  locally  too  bent  is 
reasonable,  so  that  we  then  would  use  the  Resilience  Method  to  adjust  in  such  a  way  that 
the  reasonable  bent  is  retained. 

We  adjust  locally  by  Resilience  Method  in  the  interval  once,  where 

h  *  1  or  2  for  case  1  or  2  respectively. 

C.  Adjustment  of  local  part  of  spline  which  has  one  directional  large  curvature  jump. 

This  occurs  when  c ,  c ^ ,  Cj+^,  cj+2  have  the  same  sigature  and 


>  K|cj+1l,  Cj  >  Kc^c^.  c.cj+2 


>  KCj+1 


(38) 


K  as  in  (29).  Take 


6y.  =  ao,  6y  .  -  6 a  (39) 

J  3+1 

in  (14)  and  in  general  let  a  •  2,  B  *  -1,  but  when  (Xj,yj)  Is  a  fixed  point  let 
o*0,  similar  for  S.  Then  determine  o  by  the  condition 


and  adjust. 


Kc 


1+1 
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Interchanging  and  Cj+1,  c^_1  and  cj+2*  we  obtain  another  case  needing 

adjustment. 

After  local  adjustment,  the  three  qualitative  requirements  may  still  not  be  satisfied, 
so  that  we  must  do  qualitative  adjustment  again. 

5.  A  few  words |on[  programming. 

The  effective  program  for  data  smoothing  ought  to  have  many  capabilities  such  as: 

1.  Can  rotate  axis,  such  that  the  condition  |y'(x)|  <  1  on  the  whole  spline  is 
satisfied. 

2.  Can  deal  with  nearby  points  or  not,  as  wanted. 

3.  Can  do  any  of  the  five  methods  of  smoothing  separately:  the  Resilience  Method  R,  the 
Rule  Modelling  Method  M,  the  method  A,  B,  C  of  4. 

Can  do  the  above  methods  in  some  suitable  order  such  as  RM,  MAM,  MBM,  MCM,  MABCM, 
RMABCM. 

It  is  also  required  that  using  the  program  we  can  obtain  the  smoothing  result 
quickly.  This  requires  that  the  coefficients  of  the  spline  be  calculated  quickly,  e.g.,  by 
some  local  method.  For  this,  we  derived  a  series  of  formulas,  they  are  written  down  in 
[3]. 


ACKNOWLEDGEMENT 

The  author  thanks  Professor  Carl  de  Boor  for  his  valuable  help  and  careful  reading  of 


this  paper 


REFERENCES 


[1]  (3.  Birkhoff  and  C.  de  Boor,  Piecewise  Polynomial  Interpolation  and  Approximation, 
"Approximation  of  Functions",  H.  L.  Garabedian  ed.,  Warren,  Michigan,  1964,  pp  164 
190. 

[2]  G.  C.  Dong,  D.  C.  Wang  and  X.  K.  Lin,  The  Resilience  Method  for  Mathematical  Ship 
Lofting,  Academic  Press,  1978,  Peking. 

[3]  G.  C.  Dong,  D.C.  Vang  and  X.  K.  Lin,  The  Resilience  Method  for  Mathematical  Ship 
Lofting,  Technical  Report,  1978,  Zhejiang  University. 


GCD/jv» 


a: K  •  rtf.-?-'- 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whm  Omit  Kntmnd) 


REPORT  DOCUMENTATION  PAGE 


I.  Sc  POST  NUMBER 

#2151  v 

*.  TATLE  (mnd  Submit) 

|  pATA  SMOOTHING . 


jip.C  READ  INSTRUCTIONS 

MUt- _ BEFORE  COMPLETING  FORM 

2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


1.  TWrem^Wmwruoo  covered 

^ummary/RepirU^/no  specific 
re  porfa»wj~p€nod _ 

6.  PERFORMING  ORG.  REPORT  NUMBER 


IT.  AUTHOR! •> 


/  w  _ __ 

Guang-Chang jponq 


».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Mathematics  Research  Center /University  of 

610  Walnut  Street  Wisconsin 

Madison,  Wisconsin  53706 _ 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS  ,^'7N 

U.  S.  Army  Research  Office  [// . 

P.O.  Box  12211  '  ; 

Research  Triangle  Park,  North  Carolina  27709 

U.  MONITORING  AGENCY  NAME  A  AODRESSfi/  dUimront  from  Controlling  Ollico) 


B.^jCQIji TRACT  OR  GRANT  NUMBER(AJ 

(?S‘  DAAG29-80-C-004i/' 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Work  Unit  Number  3  - 
Numerical  Analysis 

\2~  WCPORT-MATE  ~  "" 


13.  NUMBER  OF  PAGES 

17 

IS.  SECURITY  CLASS,  (ol  thtm  import) 

UNCLASSIFIED 

IS*.  DECLASSIFICATION/ DOWN  GRADING 

SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  thll  Rtporl)  “  * 

Approved  for  public  release;  distribution  unlimited. 


[  17.  DISTRIBUTION  STATEMENT  (ol  II »•  obatrmct  mntmrod  In  Block  20 ,  II  dllloroni  Irom  Rmport) 


pTi.  SUPPLEMENTARY  NOTE* 


119.  KEY  WOROS  fConiinu*  on  toworom  oido  II  nocoooory  on d  Identity  by  block  ntmbor) 


smoothing ,  spline ,  univariate,  ship  lofting 


IQ.  ^ytiUCT  (Continue  on  rovoroo  cldo  II  nocoooory  ond  Identity  by  block  nombor) 

A  new  spline  smoothing  method  is  outlined.  The  method  is  modelled  after 
the  smoothing  techniques  used  by  draftsmen  in  ship  lofting.  Full  details  can 
be  found  in  [2],  13]. 


DO  ,XTn  1473 


EDITION  OF  I  NOV  SS  IS  OBSOLETE 


«7  * 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (**««  DM* 


